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1 Abstract 

A number of experiments for measuring anisotropies of the Cosmic Microwave Background use 
7-H ■ scanning strategies in which temperature fluctuations are measured along circular scans on the sky. 

It is possible, from a large number of such intersecting circular scans, to build two-dimensional 
. sky maps for subsequent analysis. However, since instrumental effects — especially the excess low- 

| frequency 1// noise — project onto such two-dimensional maps in a non-trivial way we discuss 

. the analysis approach which focuses on information contained in the individual circular scans. 

7-H ' This natural way of looking at CMB data from experiments scanning on the circles combines 

t • the advantages of elegant simplicity of Fourier series for the computation of statistics useful for 

constraining cosmological scenarios, and superior efficiency in analysing and quantifying most of 
. the crucial instrumental effects. 

2 ■ 1 Introduction 

The exploration of Cosmic Microwave Background (CMB) is undergoing a boom in both theoretical 
■ and experimental directions. The planning of two ambitious space missions, MAP and PLANCK, to 
be launched in the beginning of the XXI-st century, is one of the sources of stimulus for research in 



this field. In addition, several long duration balloon-borne experiments are scheduled for operation 
within the next few years time. 

Many of the future CMB anisotropy experiments will perform circular scans on the sky. MAP[|, 
PLANCK0 (Bersanelli et al., 1996), and the balloon-borne experiments TopHat^j and BEAST (Lubin, 
private communication), will collect data from a large number of intersecting circles, which will then be 
merged into two-dimensional sky maps. Smaller ground-based experiments as DIABOLO (Benoit et 
al., in preparation) may scan on only a few circles. In winter 1997 for instance, a one-ring observation 
of anisotropies was done at the POM2 telescope in the french alps with the DIABOLO instrument. 

A common problem in the analysis of the data obtained by performing one dimensional sky scans is 
that all instrumental effects (other than smoothing by the beam) occur in the monotonic time domain 
of the data stream, but need to be projected onto the often complicated geometry of directions of 
observation, which is essential for the interpretation of the results. Excess low-frequency (typically 
~ 1//) noise, for instance, is an instrumental effect that is not easily analysed when projected onto 
two-dimensional maps (Delabrouille, 1997). 

In the following, we explore the relation between the Fourier spectrum of temperature fluctuations 
measured on a circle and spherical harmonic expansion coefficients on the sphere, with associated 
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uncertainties in the framework of gaussian statistics, and investigate how circular scans permit to 
analyse naturally both the signal of cosmological origin and some of the effects coming from the 
instrument. 



2 CMB anisotropies on circular scans 
2.1 Fourier spectra 

The usual expansion in spherical harmonics of the statistically isotropic CMB temperature fluctuations 
on the sky observed with a symmetric beam reads: 

oo I 

T (™)=Y1 a imBiY em (n). (1) 

g=l m=-£ 

Here, the coefficients ag m are assumed to be zero-mean Gaussian deviates with variances given by 
(|fl£m| 2 ) = Cg (on assumption of statistical isotropy of the CMB temperature perturbations, here ex- 
pressed via the ensemble averaging), Bg is the beam response function (we have assumed for simplicity 
that the beam is symmetric), and a unit vector n defines the position on the sphere. The functions 
Yg m (n) = Yg m (9,4>) are the orthonormal spherical harmonics, defined, for — £ < m < £, as 

Yg m {6A) = Vr(0)e im< t>, (2) 

with 




Vg m {0) = \ 2i A +li l , ""I! JWcosfl), form>0, 
(I + m)\ 

'Pe\m\( G ), for m < 0, 

where PJ™ are the associated Legendre polynomials. This definition of the Yg m , together with the fact 
that the temperature is real, implies ag m = (—l) m a^_ m . 

If only one ring of angular radius on the sky is being scanned, it is convenient to use coordinates 
such that the observed circle is the set of points of the sphere at constant colatitude G. The CMB 
temperature on this circle, T(G, 0), can be decomposed uniquely in the form of a Fourier series: 

- / #T(G,0)e- im ^ = ag m BgVg m {®). (3) 
!tt Jo „ , , 
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If ag m are uncorrelated gaussian random variables, as expected in an inflationary scenario (Bardeen 
et al., 1986), the same is the case also for the ring mode amplitudes a m . The corresponding Fourier 
spectrum, whose components we shall denote as T m , is the one dimensional analog of the Cg spectrum 
on the sky. The anisotropy Fourier spectrum T m is obtained from the Cg spectrum by: 

{a m a* m ,) = r m 5 mm >, (4) 

and 

oo 

r m = ]T CgBjV 2 lm {Q). (5) 
e=\m\ 

This last equation allows a straightforward computation of the ring anisotropy power spectrum T m 
given the instrumental specifications (Q and Bg) and the cosmological model (Cg spectrum as a 
function of the relevant parameters). 

Figure |l] illustrates some properties of the ring spectrum r m as compared to the usual full sky 
spectrum Cg. Observational configurations chosen for the plot correspond to 1) a single ring scan of the 



PLANCK satellite observing with the High Frequency Instrument (the large ring), 2) a single ring of 
the TopHat balloon experiment (the intermediate ring), and 3) a hypothetical small ring and a narrow 
beam configuration. One should note that the ring anisotropy power spectra preserve the dependence 
on cosmological parameters observed by the curves, and, if the resolution of the instrument is 
good enough and the ring size sufficiently large, display clearly the array of adiabatic peaks. This 
promotes the strategy of observing the CMB anisotropy on the rings to a very interesting status. 
Clearly one should be able to conduct a well guided analysis of such a data with the aid of a simple 
statistic, T m , which, due to the natural geometric reduction from the whole sky to the circle, captures 
all those features of theoretical CMB anisotropy predictions — direct reflection of physical processes 
that perturb CMB temperature in the shape of power spectrum, and the traceable dependence on 
cosmological parameters — which created high level of expectations for the full sky CMB mapping 
missions. 




Figure 1: Left panel: Full sky power spectra Ct of CMB anisotropy, smoothed with a gaussian 
beam of FWHM = 7 arcmin, in various CDM cosmological models normalized roughly to match the 
CMB anisotropy detected by COBE-DMK. All models obey the primordial nucleosynthesis constraint 
Qbh 2 = 0.015. Solid lines (top to bottom) correspond to flat, matter dominated models with Hubble 
constant h = 0.5, 0.6, 0.7, 0.8. The dashed line shows a cosmological constant dominated model with 
A = 0.7 and h = 0.8, and the dot-dashed line shows an open model with = 0.3 and h = 0.65. Right 
panel: One-dimensional power spectra of CMB anisotropy on circular scans computed for the same 
models as shown in the left panel. Three groups of curves as viewed from left to right correspond to 
different specifications (by the angular radius of the circle, O, and the FWHM of the beam) of the 
geometry of observations: the left group corresponds to a small ring of O = 4°, FWHM = 5'; the 
middle group (amplitudes multiplied by 2 to avoid overlap with the other curves) corresponds to an 
intermediate size ring of O = 12°, FWHM = 20'; finally, the right group corresponds to a large ring 
of O = 80°, FWHM = 7'. Both the left and right panels show the CMB anisotropy rms amplitude 
per logarithmic increment of the relevant index t or m, taking into account the number of degrees of 
freedom per mode (factors (2£ + 1)£, and 2m, respectively). 

Ring power spectrum coefficients T m can be viewed as the estimators of the spectrum Cg integrated 
over i with a certain m-dependent window function whose coefficients are simply: 

wt ) =V im 2 (S)B i 2 . (6) 

A number of such window functions W^ m \ corresponding to the ring configurations used in Fig. [l], are 
displayed in Fig. |2[ An important, if simple, property of these window functions, and, hence, the T m 
coefficients, is that the power at mode m is generated by only those components of sky anisotropy for 
which £ > m. This is rather important for the suborbital attempts to measure the CMB anisotropy, 



which have to cope with the atmospheric effects. It has been the case for a long time that due 
to the variation from zenith to horizon of the length of path through the Earth's atmosphere, and 
the associated systematic effects in the measurements of CMB temperature, the preferred scanning 
strategy in CMB experiments was usually along the lines of constant elevation. Can one, in a suborbital 
experiment, afford to scan on any other circle than that at constant elevation? If the atmospheric 
temperature gradient is predominantly large-scale compared to the size of the observing ring, its 
contaminating effects will be confined to low-m modes of the ring anisotropy spectrum, and high-m 
modes will be algebraically decoupled from such a source of spurious anisotropy. 
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Fi gure 2: Window functions for computation of the CMB anisotropy ring power spectrum r m , see 
Eqs. (|5|), and (0), normalised to unity at their maxima. 

To what extent would a measurement of CMB anisotropy on the ring allow us to estimate the 
parameters of cosmological model used to describe the data? To be able to answer this question we 
need first to clarify the statistical properties of circular modes of CMB anisotropy, and the properties 
of instrumental noise as projected on the ring geometry of observations. We shall see momentarily 
that the estimation of uncertainties in measurements of a m due to white noise, excess low-frequency 
noise, or instrumental systematic effects is simpler and more secure in the case of ring anisotropy than 
in the case of two dimensional maps used in C(_ estimation. 

2.2 Statistics of Fourier transform of CMB anisotropy measurements on the rings 

Here we address the issue of cosmic variance of the ring modes of the CMB anisotropy. Fourier 
coefficients of temperature fluctuations measured on the circle, a m -s, are linear combinations of a£ m -s, 
whose statistical properties they hence inherit. Specifically, Gaussian distribution of a^ m -s renders 



Gaussian distribution of a m -s. Statistical isotropy of the model CMB anisotropy field results in 
a very special property of the m-independence of variance — (|a&n| 2 ) = Ce, which makes the Ci 
spectral coefficients xhe+i) distributed. In the case of Fourier coefficients of the ring data the number 
of degrees of freedom associated with a specific m value is reduced to only two — ±m waves have 
identical variance. Hence, the cosmic variance on the T m coefficients is described by the x| probability 
distribution, independent of the value of m. 
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Figure 3: Statistics of the rms power of the CMB anisotropy measured on the ring. Gray band covers 
the 68% confidence range for each rms power coefficient and scales identically at each value of m in 
terms of the variance, which is shown by solid line. Dashed line shows the loci of the most probable 
values of the rms power, illustrating the skewness of the distribution. 

In Figure |3] we illustrate the statistical properties of ring Fourier coefficients of CMB anisotropy. 
The plot shows a few features of probability distribution of the rms power in mode m. It is clear 
that the cosmic variance in the ring Fourier coefficients is large, and, unlike in the case of full sky 
measurements of CMB anisotropy, when available at small enough angular scales, the statistics of 
power at high m do not improve compared to those in low-m modes, as the ring configuration allows 
just two degrees of freedom per mode, independent of m. It is clear that the only opportunity for 
improvement of statistics of the measurements of ring modes of CMB anisotropy rests with repeating 
the observations on many rings. The essential issues of assessing the ring to ring covariance of the 
Fourier modes of CMB anisotropy will be addressed in detail in a separate paper (Gorski, Delabrouille, 
and Hivon, in preparation). 

3 Noise on circular scans 

Now we turn to the estimation of the projection of instrumental noise onto a ring of data. Noise, 
unlike the CMB temperature anisotropies, is a time-dependent process which is independent (to first 
order) of the position on the sky. The same is true of many instrumental effects in the measured signal 
(system impulse response, sampling,...), but we first limit our discussion to a simple Gaussian noise 



as an illustration of the connection between the measured data stream and the signal reprojected on 
the ring. 

Let us assume that data is collected by a total-power experiment as PLANCK, TopHat, BEAST 
or DIABOLO, and that a circle on the sky is scanned continuously N times, with a spinning frequency 
/spin (spinning period T sp i n ). The output of each detector is a signal s(t), whose value as a function 
of time can be written as: 



s(t)=u(0(t))+n(t). (7) 

In this equation, <p(t) = 2irf SV i n t is the longitude angle on the circle. The way it depends upon time 
reflects the scanning strategy. u{(p) = T(0, 0) is the useful signal from the sky, i. e. the temperature 
anisotropy smoothed with the beam. It is a 27r-periodic function of eft for fixed 0. If anisotropies 
are correctly represented by a gaussian random field, the Fourier components u m of the Fourier series 
decomposition of u(4>) are realisations of the random variables a m defined in the previous section. 
The noise n{t) is assumed to be a realisation of a Gaussian stationary random process with zero mean 
and a bilateral spectrum S n {f) which, over a large range of frequencies, can be typically of the form: 

S n (f) = a (l (8) 
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This is the simple case of the so-called 1// noise. 

Given the signal timeline as defined in Eq. (|^), the estimation of the Fourier components of 
anisotropy on the ring can be done by computing the integral: 

1 /"NTspin / 2inmt\ 
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/ s(t)exp - — dt 
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= u m + — — / n(t) exp -— dt. (9) 

^Vispin JO \ ^spin J 

Under the present hypotheses (no foregrounds, uncorrelated gaussian noise), this quantity s m is an 
unbiased estimator of a m , which variance originates in only two sources of uncertainty. First, there is 
the theoretical "cosmic" or "sample" variance due to the fact that u m is just one realisation of a m , as 
discussed in paragraph 2.2. Second, there is the uncertainty due to the detection process (i. e. due to 
the noise n(t)). Let us denote as n m the quantity 

n m = —— / n(t) exp -— dt. (10) 

i S pin JO \ J spin / 

Each of the n m -s can be viewed as a realisation of a gaussian random variable. They obey the 
correlation statistics: 

/oo 
S n (f)H m (f)H* m ,(f),df (11) 
-oo 

where H m (f) is given by: 

H m (f) = e (*"C///M»-"0) x ^(^(///spin-m)) 

Interestingly, but not surprisingly, there may be some correlations between Fourier amplitudes of 
the noise reprojected on the ring, depending on the nature of the noise spectrum S n {f). However, if 
the number iV of scans on a circle is large and the noise spectrum S n is not too steep (1// is fine), 
the expectation value of n m n* m becomes: 

{n m n m ,} ~ — — 8 mm >. (13) 

*■ " -^spin 



Thus, for each individual m value, the uncertainty in s m Sm induced by noise depends on / sp i n 
through the noise spectrum value 5 n (m/ sp i n ), and upon the total integration time iVT sp i n . For a given 
experiment, the spinning frequency and total integration time should be optimised by comparing the 
level of the noise spectrum to the expected spectrum of the signal of astrophysical origin. Usually, the 
noise spectrum decreases with increasing /, and thus the minimisation of n m n* m requires to spin fast, 
up to a compromise which depends on other experimental considerations. 

If the number of scans on a circle, N, is large, and the noise spectrum S n is not too steep, the 
expectation value of s m s^ averaged over sky realisations of the CMB and noise realisations is: 

( Sm ,s* m ) ~ (r m + 5n ^ in) ) <W (14) 

Thus a direct comparison of the signal and noise spectra for typical next-generation experiments 
is quite straightforward in the ring formalism and is presented here in Fig. [|. 

In the top panel of Figure we have plotted together the rms values of T m for three different 
experiments and two different cosmological models, and the noise-induced sensitivity limits to indi- 
vidual T m components. Experimental parameters chosen for the plot are representative of the next 
generation balloon-borne or satellite experiments. The pure white noise contribution, proportional 
to m 1 / 2 , has the shape of an exponential because of the logarithmic scale for m. 1/f noise alone 
would appear on the plot as an horizontal line, with amplitude equal to that of the white noise at 
iTT' = "T-knce = /knee //spin- As long as mk nee is smaller than a few, the excess 1/f noise induces a slight 
flattening of the noise curve for low m, with negligible total additional noise power. 



4 Other systematics and instrumental effects 

Apart from noise, there are a few effects of instrumental origin which can be understood much easier 
on timelines than on maps of the sky. Some of these effects can be described as those of filters with 
a known impulse response, hit). The Fourier transform of the impulse response, H(f), is known in 
signal processing as the transfer function of the filter. For such filters, the filtering theorem states 
that if a signal u(t) has a spectrum S u (t), the corresponding filtered signal, «f(£) = u(t) * h(t), has 
spectrum S u (t)\H(f)\ 2 . Two important instrumental effects can be addressed with this formalism. 



4.1 Time constant of bolometers 

A typical bolometer for CMB anisotropy observations is a temperature-sensitive resistor heated by 
incoming radiation and cooled by a heat-conducting connection to a sub-kelvin thermal bath. The 
response of such a system to an impulse of incident power is: 

h(t) = - exp(-i/r) (15) 
r 

The constant r is known as the time constant of the bolometer. It is set by physical characteristics 
of the bolometer (r = C/G, where C is the heat capacity of the bolometer and G the effective 
conductivity of the link to the thermal bath). However, other bolometer characteristics, and in 
particular the sensitivity of the bolometer, depend on the same parameters. Thus, it is crucial, 
for bolometer optimisation, to understand the effect of the time constant of the bolometers on the 
spectrum of the useful astrophysical (and especially cosmological) signal. 



The transfer function H (/) corresponding to the impulse response of equation 15 is: 



(16) 
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Figure 4: Top panel: Power spectra T m of signal and noise for the following set-ups of CMB anisotropy 
measurements: 1) the ring radius of = 4° and FWHM = 5 arcmin — light lines, 2) G = 12° and 
FWHM = 20 arcmin — medium-width lines, and 3) G = 80° and FWHM = 7 arcmin — heavy lines. 
Solid lines show the predictions of standard CDM model, and dot-dashed lines those of an Qq = 0.3 
open model, both normalized to COBE-DMK. Three gray bands show the range of rms values of 
instrumental noise contributions. Bottom edges of the gray bands correspond to pure white noise 
behaviour of the detector. Solid lines inside the bands illustrate the noise enhancement due to the 1/ / 
component when the knee frequency is equal to the spin frequency on the ring scan. The upper edges 
of the gray bands correspond to /knee/Zspin = 5 (excess low frequency noise grows with /knee / /spin )• 
The assumed (white noise) detector performance is characterized by sensitivity of 1 mK^/i and total 
integration time of 20 and 100 hours for the upper two bands (light and medium gray respectively), 
while the third, narrow, dark gray noise band corresponds to a sensitivity of 20/uK-^/s and integration 
time of 2 hrs. Bottom panel: Illustration of the effects of sampling and bolometer time constant (see 
text). Heavy lines correspond to spin period of 1 min, and illustrate the effect of sampling at 10 ms 
- solid line, and bolometer time constant of 6 ms — dashed line. Medium-width and light curves 
correspond to spin period of 10 s, and 3 s, respectively. For these two spin rates the solid lines illustrate 
the effect of sampling at 10 ms, and the dashed and dot-dashed lines show the effect of bolometer time 
constant of 6 ms, and 40 ms, respectively. 



Using the relation between m and frequencies of the signal, m = /// sp i n , the attenuation function A m 
on the T m corresponding to the effect of the time constant of the bolometers is: 

A m = \H m \ = - ■ — - yz (17) 
and instead of measuring the quantities T m , the instrument measures A m T m . 



4.2 Effect of the sampling frequency 

Another instrumental parameter which needs to be set carefully by designers of CMB experiments is 
the sampling frequency. It should be high enough in order to avoid aliasing, but making it unnecessarily 
high results in a very large set of data, with the associated problems of data storage and telemetry, 
especially for satellite or long duration ballooning experiments. 

Again, the effect of the sampling rate can be understood very simply as that of a filter. If sampling 
is modelled by a perfect integrator over a period To (which is the inverse of the sampling rate) for 
obtaining each sample, the corresponding impulse response is: 

h(t) = ^rect(To), (18) 

where rect(To) is the function whose value is 1 between t — To/2 and t + To/2, and elsewhere. The 
corresponding transfer function is 

= sj^pi 

and the corresponding attenuation function A m on the T m is: 

i l9 sin 2 (7rmf SD ; n To) , . 

Am = \H m \ 2 = , V , ^2 20 

Attenuation functions for these two effects are plotted in the bottom panel of Figure [|. 



5 Conclusion 

In this paper, we presented the formalism to project the CMB anisotropy predictions and instumental 
noise effects on the rings of data. Because of the feasibility of experimental setup for CMB observations 
with circular scans, and simplicity of the data analysis on such scans, we believe that collecting 
CMB anisotropy data in this format is an interesting option for the next generation of anisotropy 
experiments. It provides a natural frame for studying the statistics of CMB anisotropies, in which 
many instrumental effects can be modelled and analysed under the most natural connection between 
the time stream of data and the spatial distribution of the directions of measurements. We attempted 
to illustrate this point comprehensively in our Figure ||. Given that this comparative plot allows to 
present simultaneously a natural rendition of a number of factors that affect the attempts to measure 
the CMB anisotropy on the celestial rings, we believe that Fig. |] should provide significant insights 
on optimization of experimental setup for future CMB experiments. 

For the most ambitious CMB programmes such as PLANCK, many such scans will have to be 
merged into two-dimensional maps. The circular scan method described in this paper is a useful 
intermediate step in the analysis of such data in the global analysis of CMB anisotropies on the 
celestial sphere. 
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